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SUMMARY 

Methods are presented that use general correlative time- 
response Input and output data for a linear system to determine 
the frequency- response ftmctlon of that system. These methods 
give an exact description of any llneeo? system for which such 
transient data are available. 

Examples are shown of application of a method to both an 
underdamped and a critically danq>ed exact second-order system, and 
to 6tn exact first-order system with and without dead time. Esqperl- 
mental data for a turbine-propeller engine showing the response of 
engine speed to change In propeller-blade angle are presented and 
analyzed. 


INTRODUCTION 

A basic problem confronting the designer of engine controls 
Is that of determining the behavior of the complete engine and 
control system xmder veirylng conditions of operation. This ques- 
tion becomes peurtlculeirly acute when veirlous engine parameters, 
for example, shaft torque or turbine- Inlet temperature, must be 
closely controlled In order to prevent engine damage or even fail- 
\ire. The problem then becomes one of matching transient behavior 
of the control to that of the engine In' order to attain a desired 
system response. 

General methods do not exist for the solution of the equa- 
tions of motion of nonlinear systems. As a result, the analysis 
of such systems may be ln 5 >ractlcably difficult. If, Indeed, a 
solution can be found at all. Because the behavior of many non- 
linear systems may be satisfactorily approximated by the assiunp- 
tlon of system linearity and because general mathematical methods 
and techniques for the analysis of linear systems are readily 
available and (in comparison with present nonlinear methods) rela- 
tively simple, to apply, the asstmq>tlon is generally that the 
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system being studied Is llneeu:. The methods of this report are 
based on such an assumption. 

Approaches utilized In the analysis of linear systems (refer- 
ence 1, pp. 17-18) are: (l) transient analysis, In which the char- 

acteristic time response of the system Is determined for standard 
Inputs such as step and Impulse functions, or (2) frequency anal- 
ysis, In which system behavior Is characterized by the steady-state 
amplitude and phase relations of the system Input and output for 
sinusoidal Inputs of various frequencies. 

For Inputs such as step or Impulse functions, the Inherent 
system characteristics might be obtained by fitting an equation to 
the output function and from it deriving the differential equation 
of the system. If the Input and output functions were of any gen- 
eral form, fitting differential equations to the data might still 
be possible, but for systems of inherently high order the accuracy 
of such a procedxure would be low. In general, use of the time- 
response form of system description involves dealing with convolu- 
tion Integrals (reference 2, p. 54) when the xmlt so described Is 
to be combined with other units. For complex systems, this descrip- 
tive form Is not readily adapted to manipulation. 

A linear system is known to be characterized by Its steady- 
state response to all frequencies of sinusoidal inputs. The 
frequency-response form Is useful for the description of linear 
systems because of the ease with which various system character- 
istics can be manipulated In dealing with combinations of units. 

For example, the over-all anplltude ratio of output to Input for 
any frequency of input to a system consisting of several tmits in 
series Is obtained by a simple multiplication of the amplitude 
ratios for the several units. A possible diffictilty in the applica- 
tion of this method to analysis of engines is the necessity for 
maintaining a sinusoidal Input of constant amplitude euid frequency 
for a length of time sufficient to insure disappearance of transient 
effects in the output. 

The nature of real physical systems may make actual Imposi- 
tion or measurement of step or Impulse Inputs Impossible; In addi- 
tion, sinusoidal inputs may prove Inpracticable . The need thus 
Gurlses for feasible methods capable of handling data of any general 
form. Such methods that xise general correlative time-response data 
for system Input and output to determine the frequency-response 
characteristics of the system have been developed at the KACA Lewis 
laboratory and are presented in this report. These methods give an 
exact description of emy linear system for which such data are 
available. 
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Three exact methods of obtaining the frequency-response func- 
tion of a system are shovn for system-equlllbrlvan final conditions* 
A graphical approximation to one of the methods Is also glran. 
tfodlflcatlon of one of the exact methods to accoxmt for oscilla- 
tory final conditions is shown and treatment of dead time Is pre- 
sented. As an illustration of one of the methods^ three exanqtles 
are given, two based on analytically determined system-response 
curves and one on experimental data obtained at this laboratory 
for a turbine-propeller engine. 
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SYMBOI^ 

The following symbols are used In this report: 
amplitude of system steady oscillation 
constants 

system frequency-response function 

system transfer function 

general function of time 

first time derivative of f(t) 

system transfer function including dead time 

Laplace transform 

complex number 

anq>lltude ratio of F(iO)) 

time, seconds 

general time- dependent Input 
general time-dependent output 
n^^ time derivative of y 

frequency of system steady oscillation, radians per second 

/ 

system dead time, seconds 
damping ratio 
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6 phase angle of F(1 (d), radians 

T system time constant, seconds 

^ phase angle of system steady oscillation, radians 

CD frequency of system sinusoidal Input, radians per second 

Subscripts: 

! 

0 last texm of differential equation 

f final value 

k general tezm of series 

m degree of differential eqiiatlon of Input 

n degree of differential equation of output 

r last term of series 


ANALYSIS 

The system considered In the follovlng derivation Is assumed 
to be linear. In particular, behavior of the system Is assumed 
representable by a linear differential equation with constant 
coefficients . The block diagram of a system having an Input x 
and a corresponding output y Is shown In figure 1, ^ The sym- 
bol F(p) Inside the box represents the system operator which, 
acting on the Input x, yields the output y. The linear differ- 
ential equation relating x and y may be written as 


®n ^ + ®n-l 
dt’^ 


dt^-1 




dt“^ dt®“^ 


. + bQX 


( 1 ) 


If the system Is Initially at rest (or In equilibrium), term- 
wise application to equation (1) of the Laplace transfoz^tlon 
(reference 2, pp. 51 to 56), defined as 

nco 


Ljf(t)] 


f(t)e dt 


( 2 ) 
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and factoring of the result give 


(®nP“ + + • • • + ao) I*(y) = + 


m-1 


or 




n xx-1 

»nP + ^n-lP + 


Equation (3) may be formally written as 

L(y) a F(p) L(x) 


®0 


where 


+ bQ)L(x) 

(3) 


( 4 ) 


F(P) 


^inP ^^-iP + • • . + bQ 

v“ Vi^“"^ + . . . + sq 


The transfer function of the system F(p) is defined as the ratio 
of the Laplace transform of any normal, response of the system to 
the Laplace transform of the input producing that response. Normal 
response is the response of the system when initially at rest 
(reference 2, ,p. 26), 


From equations (2) and (4), the system transfer function may 
be written as 



Equillbrltmi Final Conditions 

Derivation of frequency-response function. - The f requency- 
response function F(iW) is formally obtained directly from equa- 
tion (5) by replacing p by iO) (reference 1, pp, 96-98). Then, 
for a sinusoidal input x of frequency (O, the output y 
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ultimately is sinusoidal at the same frequency hut vith a relative 
magnitude and phase angle equal to the magnitude and the phase 
angle of the coiBplez number F(i(i)). The frequency- response func- 
tion is defined as 



( 6 ) 


The term 


-i(Ct 

e 


is oscillatory and the Integrals of equation (6) 


may not converge unless 




converge. 


These conditions impose a severe restriction on the types of func- 
tion that may be used in equation (6) because of the lnq>lled 
requirement that z and y vanish as t Increases vlthout limit. 


As vill be shovn, the previous restriction on z and y may 
be removed by suitable modification of equation (6). A much vlder 
selection of input and output functions is thus permitted. 


For a system initially at rest, 

L[f(t)] »pL[f(t)] 

An alternate erpression for the transfer function then is 
or, from eqiiatlons (2) and (8), 



(7) 

( 8 ) 


The frequency-response function then becomes 



NACA TN 1935 


7 


F(10)) 



(9) 


when p is replaced by l(i). 

For any input or output function ending in equilibrium, x' 
and y* become and remain zero for t equal to or greater than 
some time t^. The integrals of equation (9) therefore converge. 

From Euler's relation. 


-icot 

e 


cos (Ot 


i sin (Ot 


equation (9) may be written 


F(i(o) 




sin (i)t dt 


sin (Ot dt 


( 10 ) 


This form of the frequency-response function can be used with any 
data for which the input and the output begin and end in a steady 
state. 

The integrals of equation (lO) can be evaluated in various 
ways. One theoretlceilly exact method, which has been used at the 
HACA Lewis laboratory, utilizes a rolling- sphere harmonic analyzer. 
This type of analyzer, the same in operating principle as the 
device described in reference 3, was designed to obtain the Fourier 
coefficients from cyclic data, but the manner in which the coeffi- 
cients are determined results in evaluation of Integrals of exactly 
the form of those in equation (10). The frequency- response func- 
tion can thus be obtained for any frequency by operating the 
analyzer over the output cxurve to evaluate the integrals in the 
numerator and over the input cui*ve to determine the integrals in 
the denominator of equation (lO) without replotting the data in 
any other form. 

The frequency spectra for nonperiodic phenomena are determined 
in reference 4 by evaluating Fourier Integrals with a rolling- sphere 
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hanoonlc analyzer In a manner similar to that followed in this 
report in connection with equation (10) . The method of reference 4, 
however, differs from that of this report in being restricted to 
functions that begin and end at zero. In addition, the concepts 
of a system transfer or a system frequency-response function are 
untreated. 

Modifications of frequency-response derivation. - The frequency 
response may, at times, be desired directly in teims of the input 
and output ftinctlohs themselves rather than their time derivatives. 
For example, if an electronic device were used for analysis, the 
avoidance of the added conplexlty and attendant inaccuracy of dif- 
ferentiating circuits ml^t be desirable. 

If both X and y come to a steady state by or before some 
time tf, equation (9) may be rewritten as 



because for t > t^ the Integrands of equation (11) are zero. 
Integrating equation (11) by parts gives 


F(iO)) 


... I -i(0t -i(i)t 

iW / ye dt + ye 

Jo 


10 ) 


-i0)t -iwt 

xe dt -I- xe 


Jo . 


tf 

0 


or 


F(i0)) 


^ sin 0)t^\ / / 

y cos (Ot dt - y^ — j - i^ y sin 0)t dt - y^ ^ 


cos cot^i 


X cos (Ot dt - Xf 


sin cot^' 


(0 


- l{ y X sin cot dt - Xf 


cos (Ot^ 


(0 


( 12 ) 
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An alternate procedure, leading to a somevhat sljiqsler ejq>res- 
slon, en^>lo7S the principle of superposition. For example, 7 may 
be considered the linear combination of a step function vhose 
ordinate is jf and a second function vhose ordinates are those 
of 7 displaced downward by the value of y^ . Such a separation 

of y into its conponent functions is shown in figure 2. An ana- 
lytical expression of the process may be derived as follows: 


The transfer function of the system is 



Becaxise of the linearity of the transformation 

ow 


F(p) = 


r 


(y-7f)e"^^ dt ■»• / yfe"^*' dt 


-pt 


J (x-Xf )e'^* dt + 


rv- 


dt 


(5) 


03 


(y-7f )e'P* dt + ^ 

f (x-Xf )e"^^ dt + ^ 


Setting p s i£i) then gives 


, , -iO)t 7f 

(y-y,). dt + ^ 

, . - icot ^ 

(x-Xf )e dt + — 


F(i(0) 
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F(i&)) 


n 00 

/ (y-Jf ) cos (i)t dt - 1 

(y-7f ) sin (Dt dt -I- ~ 

n 00 

J (z-3^) cos (Dt dt - 1 

1 ix-Xf) sin (Dt dt + ^ 


(13) 


Equation (13) defines F(iu>) for all frequencies except (o sO and 
F(0) can readily be shown to be y^/x^ . 


Approximation of frequency-response function. - If the Input 
and output functions are approximated by a series of straight- line 
segments, as shown In figure 3, then the Integrals of equation (10) 
may be evaluated by escpanslon in finite 

derived as follows: From figure 2, the 
may be approximated by 


sine and cosine series 
Integral 


y'e dt 



(14) 

If the fltst term of the first sumnatlon In equation (14) Is evalu- 
ated and the Index on the second summation Is shifted to k-1, the 
ej^resslon becomes 
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for = 0 . 

If the t's are evenly spaced, 

^ t 

*k+l - H 

Bence, 

-5LLiii 


or 


y‘k - y’k+i 


2yk " <yk-i 7k+i) 


Then 



Euler's relation, 

-i(«)t 

e s cos (Ot - i sin (ot 


applied to both members of equation (15), gives 
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from vhlch 



(17) 


where » t^, 

A Blmllar pair of espressions is obtained for the input func- 
tion, with z merely replacing y throughout. 


Oscillatory Final Conditions 

The procedures thus fea* described epply only idieh the system 
comes to rest at some new equilibrium condition. If the output 
does not reach equilibrium but continues to oscillate about some 
mean value, a modification in the method is necesscury in order to 
account for the oscillatory cosqpohent. 


System linearity permits resolution of the output y into 
two coBiponents y^^ and yg, whose algebraic sum is the original 

output, as shown in figure 4, The block diagram equivalent to this 
resolution of y into components is given in figure 5. The 
outline indicates the over-all transfer function of figure 1. ¥r<m 
figure 5, F(p) may be represented as the linear sum of two trans- 
fer functions, Fj^(p) and Fg(p), operating in paznllel. 

The following equaitlon may then be written: 

F(P) = Fi(p) + Fg(p) 
or 


F(p) 


L(yi’) L(yg') 
L(z«) 


i-Cyi') + p L(y2) 


( 18 ) 
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If the steady osclUatlon Is slnqsle harmonic, the oscillatory com- 
ponent of the output may be ea^ressed as 


yg = A sin (3t+^) 

for t > 0. Then 

p L(yg) = 

p +3 

and. 


F(p) 




dt + (3 cos 0+p sin 0) 

P^-t-3^ • 



(19) 


When p Is replaced by 1(0 In equation (19), the frequency- 
response function becomes 


F(1(0) = 



Treatirant of Dead Time 

Dead time In a linear systmn Is the time difference between 
the Initiation of a disturbance to the system and the beginning 
of a response to that disturbance. The mathematical equivalent of 
dead time then Is a trc^latlon of the output function along the 
time eo:ls In a positive direction; that Is, If 

L[y(t)] = F(p) L[x(t)] (4) 

for a system without dead time, then for the system with dead time 

L|y(t-At)J « G(p) L[x(t)j - (21) 
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vhere G(p) is different from F(p). The relation hetveen G(p) 
and F(p) is x^adily determined by an e^qpansion in a Taylor's 
series of y(t-dt) in the neic^orhood of t: 

y(t-At) B y(t) - Aty'(t) + ^ y"(t) - ^ y"*(t) + . . . 

6 e • w • 

Then 

L|>(t-At5 B L|^(tg - At L|>'(tg + ^ L^(tj - ^ L|^"'(t)]+ . . . 

For initial conditions of equilibrium 

L ^(t)] = p*^ L [y(t)] 

Hence^ 

L|y(t-At)]BL|^(t)]- AtpL[y(t)] + ^^L|y(t§-^^L(V(t)]+ . . . 

“ (i ■ Atp + + . . .^ I- [y(t)] 

« e"^^ L[y(t)] 

From equations (21) and (22) 

I*&-(t)] = e^^^ G(p) L |x(t)] 

Conq>arison of equations (4) and (23) shows that 

G(p) » F(p) 
or 

G(p) B e F(p) 

from which the corresponding frequency-response functional relation 
is 

G(l(o) B p(lco) (24) 


( 22 ) 

(23) 
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It can easily be sbovn that the effect of dead time on the 
frequency-response function Is to decrease the phase a^le alge- 
braically by the amount ttoa, leaving the amplitude ratio unchanged 
In general, F(l(i)) is a complex number and may be e^qpressed as 

P(i<D) = Ee^® 

From equation (24) 

G(1(0) = e“^^e^® 
or 

O(i(0) . 


EXAMPLES 

Examples of the application of the method of this report to 
determination of the frequency-response function of seyereO. sys- 
tems are found in figures 6 to 8. The use of the function thus 
found has not been covered because the detail vlth which such use 
has been treated elsewhere (for example, in references 1 BrnU 2) 
makes extensive treatment in this report unwarremted. 

Examples 1 and 2 Illustrate the types of frequency response 
obtained from several different types of time function. A method 
of this report was used with the specified time functions to deter- 
mine the points shown in figures 6 and 7. Frequency- response func- 
tions determined directly from the known transfer functions of the 
systems used were also plotted for comparison. Example 3 presents 
data for a turbine-propeller engine and shows the type of frequency 
response obtained, using a method of this report. 


Exaople 1 

A second-order system haying the transfer function 

F(p) = -5 (25) 

p + 2^p + 1 

is shown in figure 6 for both the tmderdamped and critically damped 
cases. The frequency-response function of the system is obtained 

directly from equation (25) by replacing p by io). For ' t = i 

2 
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F(10)) 


-(I? + 1(1) + 1 


1 - (£T 


- 1 


(0 


1 - (o2 + 1 - (D^ + (0^ 


(26) 


and for t = 1 


F(i(0) = 


-a:? + 21(1) + 1 


1 - 0)^ _ 2(1) 


(l+(i?)2 (l+d?)2 


(27) 


System response to a unit step input is given for t = 3 


;5 


(28) 


and for $ » 1 by 


1 -e"* (1 + t) 


(29) 


Frequency-response curves calculated from equations (26) and (27) 
are shown In figures 6(c) and 6(d). The rolling- sphere analyzer 
was used on the transient-response curves detemined from equa- 
tions (28) and (29) to obtain from equation (10) the points shown 
on the frequency-response cxnrves of figure 6. Because a \mlt step 
input was used, the dencnnlnator of the rl^t-hand side of equa- 
tion (lO) is equal to 1 for all values of frequency and F(i0)) can 
be obtained by operation of the analyzer over the output curves 
only. 


Example 2 

A first-order system is shown in figure 7. 
no dead time, the transfer function is 


If the system has 


1 + Tp 


F(p) = 
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and the corresponding frequency-response function Is 

F(ico) = 


1 + Tlw 


1 _ i _ 


1 + (0 ^ 1 + 0 )^ 


(31) 


For a unit step Input the output Is 


t 

y = 1 - e ^ (32) 

If the systCT has a dead time At, the transfer function hecomes 

-Atp 


G(p) = 


e 


1 + Tp 


(33) 


and 


G(1(o) = e F(Ki)) 


(34) 


The response to a unit step Input then Is 


t-At 

T 

y = 1 - e 


(35) 


for t > At. The frequency- response curves In figures 7(d) and 
7(e) were calculated from equations (31) and (34). The func- 
tion F(l(o) was found from equation (31) and G(1 Cd) was then 
determined from F( 10 )) , by Increasing the phase angle of F(l0)) 
negatively by the amount Ato), leaving |P(l(o)| unchanged. The 
rolling- sphere analyzer was vised on the output curves of equa- 
tions (32) and (35) to obtain from equation (lO) the points shown 
on the frequency- response curves of figure 7. The use of a unit 
step Input to the system makes the denominator of the right-hand 
side of equation (10) equal to 1 for all values of frequency. The 
frequency-response function can therefore be obtained by operation 
of the analyzer over the output curves only. 
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Example 3 

Experimental transient data obtained at this laboratory for 
the speed response of a turbine-propeller engine to changes in 
propeller-blade angle are shown in figure 8(a). Because the 
frequency- response function is expressed as a ratio, \mits of the 
variables involved are unin 5 )ortant and the ordinate of the data 
curves has been calibrated in relative values, the difference 
between initial and final values of the curves being set equal to 
unity. 

The frequency- response function F(iO)) of this engine for 
response of engine speed to changes in propeller-blade angle was 
found from equation (10). The rolling- sphere analyzer was used 
on the data. In figure 8(b), F(iw) is shown as a frequency- 
locus plot in the complex plane; figure 8(c) is a plot of ampli- 
tude ratio and phase angle against frequency. The closeness with 
which the data approach a semicircle in the complex plane indicates 
that the response of the system is substantially first-order. This 
inference is further borne out by the amplitude-ratio plot. The 
asynqjtotes and the solid curve were calculated for an exact first- 
order system. 


SUMMAEY OF RESULTS 

Methods are presented for the determination of the frequency- 
response function of any linear system from general correlative 
time-response data for input and output of that system. Equations 
were developed for system equilibrium and oscillatory final condi- 
tions and for the effect of system dead time on the frequency- 
response function. 

For illustrative purposes, a method of this report was applied 
to several linear systems for which the frequency- response functions 
were known. 

Experimental data obtained at this laboratory for the speed 
response of a turbine-propeller engine to change in propeller -blade 
angle are presented and analyzed. 


Lewis Fli^t Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 11, 1949. 
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X — ►- 


P(p) 


Figure 1. - Block diagram of linear system. 



(a) System-output time function. 



(b) Components of system-output time function. 

Figure 2 , - Separation of system-output time function Into 
components for equilibrium final conditions. 
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(a) System-output time function. 




(b) Components of system-output time function. 

Figure 4. - Separation of system-output time function into com- 
ponents for oscillatory final conditions. 
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Figure 5. - Block-diagram representation of 
linear system for oscillatory final con- 
ditions. 
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1 



P^ + 2tp + 1 



(a) Second -order system representation. 



t , sec 

(b) Response of system to unit step Input, 

Figure 6. - Frequency-response function for second-order system. 



Phase angle, deg Amplitude ratio Imaginary part of F(l(o) 
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Real part of P(l(i)) 



(c) Frequency-response function In complex plane. 



w, radlans/sec 



(d) Amplitude ratio and phase angle of 

frequency-response function against w. 


Figure 6. - Concluded. Frequency-response 
function for second -order system. 
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(a) Representation of exact first-order 
system with no dead time. 



(b) Representation of exact first-order 
system with dead time. 



At 1 »- 

T 0 1 2 3 4 5 6 



(c) System response to unit step Input. 

Figure 7. - Frequency-response function for first-order 

system. 




Phase angle, deg Amplitude ratio ^ Imaginary part 
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Real part of Pdco) 



A Without dead time 
V With dead time 


plane . 




(e) Araplitud# ratio and phase angle of 

frequency-response function against tco. 

Figure 7. - Concluded. Frequency-response 
function for first-order system. 








Phase angle, deg Amplitude ratio 
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(c) Amplitude ratio and phase angle of turbine- 

propeller frequency-response function against to. 

Figure 8. - Concluded, Frequency-response function of 
turbine -propeller engine. 
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